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BRANCHED MICROSTRUCTURES IN THE GINZBURG-LANDAU 
MODEL OF TYPE-I SUPERCONDUCTORS 

SERGIO CONTI, FELIX OTTO, AND SYLVIA SERFATY 


Abstract. We consider the Ginzburg-Landau energy for a type-I superconductor in the 
shape of an infinite three-dimensional slab, with two-dimensional periodicity, with an 
applied magnetic field which is uniform and perpendicular to the slab. We determine the 
optimal scaling law of the minimal energy in terms of the parameters of the problem, 
when the applied magnetic field is sufficiently small and the sample sufficiently thick. 
This optimal scaling law is proven via ansatz-free lower bounds and an explicit branching 
construction which refines further and further as one approaches the surface of the sample. 
Two different regimes appear, with different scaling exponents. In the first regime, the 
branching leads to an almost uniform magnetic field pattern on the boundary; in the 
second one the inhomogeneity survives up to the boundary. 


1. Introduction 

Superconductivity, discovered in 1911 by Kamerlingh Onnes, is a phenomenon happening 
at low temperature in certain materials which loose their resistivity and expel an applied 
magnetic held. The latter is called the Meissner effect. More precisely, when the applied 
magnetic field is small, the sample is everywhere superconducting and completely expels 
the magnetic field, while when the magnetic field becomes larger, it partially penetrates 
the sample via regions of normal phase where the material is not superconducting. If the 
magnetic field is further increased, then superconductivity is completely destroyed and the 
sample behaves like a normal conductor. 

The standard model for describing superconductivity is the Ginzburg-Landau functional, 
which was introduced in the 1950’s by Landau and Ginzburg on a phenomenological basis. 
It was later justified based on microscopic quantum mechanical principles via the Bardeen- 
Gooper-Schrieffer (BGS) theory, which explains superconductivity through the appearance 
of “Gooper pairs” of superconducting electrons. The Ginzburg-Landau model is a formal 
limit of the BGS model, and this derivation was accomplished rigorously in the recent work 
[FHSS12] . 

The Ginzburg-Landau model describes the state of the sample via a complex-valued 
order parameter u. The squared modulus of u represents the local density of the Gooper 
pairs of “superconducting electrons”. In other words, p := \u\^ indicates whether one is in 
the normal phase p ~ 0, or in the superconducting phase p ~ 1. The transition between 0 
and 1 happens within relatively thin interfacial layers (or walls). The order parameter u is 
coupled with the magnetic vector potential A, which yields the magnetic field B := V x A 
induced in the sample. The Meissner effect can be roughly understood as the fact that the 
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magnetic field B can only exist in the normal phase p = 0, or in other words 

( 1 . 1 ) pB^Q. 


Another important property of superconductors is flux quantization. If we consider a closed 
circuit well inside the superconducting region (on a scale set by the penetration length A), 
the contour integral of A has to be an integer multiple of 27r, in units of h/2e. This arises 
because Vu is very close to lAu and |m| is very close to 1. Correspondingly the flux of the 
magnetic field B through any surface with such boundary is quantized, in the sense that it 
has to be an integer multiple of 27r (again, and for the rest of this paper, in units of h/2e). 

The Ginzburg-Landau functional in a three-dimensional region Ql,t '■= (0,T)^ x (0,T) 
can be written, after appropriate non-dimensionalization, as 


( 1 . 2 ) Eqi\u,A] : — 
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Here := V — iA denotes the covariant gradient, i?ext := ^ext^a is the applied magnetic 
held, which is assumed to be uniform and vertical. The constant k > 0, usually called the 
Ginzburg-Landau parameter, is the ratio of the “penetration length” (of the magnetic held 
in the sample) A and the “coherence length” For a general presentation of superconduc¬ 
tivity and the Ginzburg-Landau model, we refer to the standard physics textbooks, such 
as |Tin96( IDG661ISJST69] . For further mathematical reference on the Ginzburg-Landau 
functional, one can see for example [SSOTj. 

Superconductors are usually classihed in type-I and type-II superconductors, according 
to whether k < 1/a/2 ot k > 1/a/2. In type-II superconductors, there is an intermediate 
regime, for low applied magnetic helds, where the penetration of the magnetic held happens 
along very thin vortex hlaments, carrying an integer hux, and around which the sample 
is normal — this is called the mixed phase. The size of the vortices is only limited by 
the hux quantization condition, and indeed in most situations each of them carries exactly 
one quantum of hux. This is the regime studied in details in |SSn7] in dimension 2 and in 
[BJQS13] in dimension 3. By contrast, in type-I superconductors the ratio of characteristic 
lengthscales k does not allow these vortex-hlaments to form and larger regions of normal 
phase appear, separated interfaces (called walls) from the superconducting phase. Each 
normal region in this case carries a magnetic hux much larger than the hux quantum. We 
will be interested only in the latter situation, and we will assume that k, is small enough, 
and also that the applied held 6ext is much smaller than the critical held, which in the 
present units is k\/2. 

The pattern of the normal phase arises from the competition of diherent ehects. On the 
one hand, the interfacial energy favours a coarse structure in the interior of the sample. On 
the other hand, the magnetic energy outside the sample favours a hne-scale mixture close to 
the interface. Therefore, the optimal pattern is expected to branch, as predicted by Landau 
back in 1938 |Lan381 ILan43] . This permits to combine a coarse pattern in the interior with 
an induced magnetic held V x A almost aligned with Bg^t at the surface, see Figure [T] for 
a sketch. Experimentally, this is manifested by complex patterns observed at the surface 
of the sample |Pro07l IPGPPOSl IPHG081 IPH09j . This phenomenon of domain branching 





























Figure 1. Sketch of the flux patterns. Left: the regime with uniform flux on 
the surface, with energy proportional to LF'^ which is the optimal 

scaling if 6ext > Right: the regime with flux concentration on 

the surface, with energy proportional to ^ which is the optimal 

scaling if 6ext < 


occurs also in other areas of materials science, as for example ferromagnetic materials, 
where the magnetization pattern is constrained to oscillate between two opposite vectors 
[Lif44l lHub67] , and in shape-memory alloys, where the strain can oscillate between hnitely 
many values, corresponding to the different martensitic variants. The average behavior 
of these branched patterns can be characterized via scaling laws: one determines how the 
minimal energy per cross-sectional area scales with the various parameters of the system, 
and shows that the optimal scaling of the energy can be achieved with branching-type 
patterns. This is usually rigorously established by showing ansatz-free lower bounds and 
complementing them with the construction of explicit branching patterns whose energy is 
estimated to have the same order in the parameters as the lower bound. This was achieved 
in martensites in |KM92l IKM94t IConOOt 100091 IC012| IZwil4| ICC 15] and in magnetic 

Te09]. 
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In the case of type-I superconductors, a similar program was carried out in |CKOn4] 
ICCKOn^ for a simplihed model: it is a “sharp-interface” version of the Ginzburg-Landau 
functional, where the order parameter u is only represented via its modulus p, which in 
turn is only allowed to take values in {0,1}; at the same time the kinetic energy is replaced 
by a constant times the BV norm of p, i.e., the perimeter of the set where p = 1, see 
Section 2.4 below for details. This resulted in a full characterization of the phase diagram 
at the level of energy scaling, and in particular led to the discovery of a new phase for very 
small applied helds, see Figure [T] for a sketch. 

We study here the full Ginzburg-Landau model, as given in (1.2), and determine the 
scaling of the minimum energy per cross-sectional area in dependence of the problem 
parameters k, 6ext and T. We prove that the energy scaling is characterized by the same 
two regimes which had been found for the sharp-interface functional. In fact, some of 
the ideas of proof of |GKO04[ IGGKO08] carry over to the full model, once an appropriate 

































































splitting of the energy, involving a “Bogomoln’yi operator” has been performed (cf. Section 
[^. The treatment of the lower bonnd contains several additional difhcnlties, mainly due 
to the fact that the Meissner condition (1.1) is only true “on average”, in the sense that 
an appropriate weak norm is small. The constructions in the upper bound, at the same 
time, need to take into account the quantization condition, locally in each tube, and to 
construct an order parameter with diffuse interfaces. The relationship with the simplified 
model is discussed in more detail in Section 12.41 

We work in an infinite periodic slab geometry, which is enough to understand the main 
surface branching features that we are interested in, hence the choice of working in the 
domain Ql^t ■= Ql x (0,T) with Ql '■= (0,L)^, with L very large, and with horizon¬ 
tally periodic boundary conditions. We recall that the Ginzburg-Landau functional is 
invariant under gauge-transformations: two configurations [u, A) and {u, A) are called 
gauge-equivalent if there exists d* G such that 


f u{x) = 

\ A{x) = A{x) + V<h(a;). 


The physical quantities are gauge-invariant (i.e., invariant under a gauge-transformation). 
This includes the magnetic held V x ^4, the energy, the density p = \u\^ and the supercon¬ 
ducting current dehned as 

(1.3) jA '■= Re(— toVam) = - [—iuVAU + iuWAu) ■ 

We will work in the space x M; C x M^) dehned to be the set of (n. A) G 

such that Eqi\u,A\ is hnite and for every fc G x {0}, («(■ -|- Lk),A{- + Lk)) is gauge- 
equivalent to {u,A) (this periodic setting was rigorously formalized in |Ode67j . see also 
[Put] ). All gauge-invariant quantities, such as p, Ja and B, are then Q^-periodic. We 
stress that periodicity is only assumed in the hrst two variables. We will also call such 
pairs {u, A) G admissible. 

Our main result, characterizing the energy in the regime of small applied helds 6ext and 
large and thick enough samples, is 


Theorem 1.1. For any K,bext, L,T > 0 such that 


bey^tL^ e 27 rZ, 


8^ext ^ ^ ^ 2 ’ 


kT > 1 


and L is sufficiently large (in the sense of (4.2)j, one has 

Egl[u, A] - («:\/26ext - bl^fL^T ~ min if 


(1.4) 


mm. 

-^-^per 


The result will follow from Theorem 3.15 and Theorem 4A below using Lemma 2A to 
separate the bulk contribution. The notation a ~ & means that a universal constant c > 0 
exists, such that c~^a < b < ca. 

The scaling result (1.4) is in the end the same as in |CKO041 ICCKO08] . after some 
rescaling of the lengths and magnetic field intensity. As in those works, the minimum 
in the right-hand side reflects the fact that two types of construction are needed, one 



















corresponding to the regime where and one corresponding 

to the opposite case. In both cases, the construction that gives the optimal scaling law is 
that of a self-similar branching tree of normal region (where the magnetic held penetrates) 
which is symmetric with respect to the X 3 = T/2 plane, and rehnes further and further as 
X 3 approaches 0 and T. The optimal “opening ratio” of the self-similar tree depends on 
the parameters of the problem and has to be chosen differently in the two regimes above. 

A hner analysis in the asymptotic regime 6 ext —t 0, including T-convergence to a reduced 
model with energy concentrated on lines, will be discussed in |CGOS] . 


This paper is organized as follows. In Section we show how the bulk contribution 
energy can be algebraically separated, via the Bogomoln’yi operator, and we dehne in 
(2.3) the functional E on which we shall focus for most of the paper. In Section]^ 


we 


prove the lower bound, hrst for the sharp-interface version of the problem, and then for 
E. In Section we prove the corresponding upper bounds; again we first work on the 
sharp-interface problem and then extend the upper bound to the full Ginzburg-Landau 
functional. 


2. Preliminaries 

2.1. Preliminary on notation. We use a prime to indicate the hrst two components 
of a vector in and identify with x {0} C Precisely, for a G we write 
a' = (oi, 02, 0) G C given two vectors a, 6 G we write briehy a' x b' = {a x 6)3 = 
(a' X 6 ') 3 . 

We shall denote sections of Ql,t by Q{z) = Ql x {z}, for integrals over Q{z) we write 
dx^ instead of dV?{x'). In the entire paper we let k, 6 ext, Li T be positive parameters which 
obey 

(2.1) 6 ext ^ ; iz < -, and kT > 1 . 
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In many parts we shall additionally require the quantization condition 6 ext-b^ G 27rZ. 

By a < 6 or 6 > a we mean that a universal constant c > 0 exists (which may change 
from line to line but does not depend on the parameters of the problem) such that a < cb. 
By a ~ 6 we mean a <b and b < a (with two different implicit constants). 

We denote by the space of traces of Q^-periodic functions u G x 

(— 00 , 0 )) with Vu G T^, and use the homogeneous norm ||m||hi/ 2 (q^) := inf{|| Vn|| 2 , 2 (Q^x(-oo,o 
where the inhmum is taken over all possible extensions. We denote by its dual 

space. 


2.2. Separating the bulk energy. Our hrst step is to subtract the bulk contribution 
to Eqi^, which will lead us to the dehnition of the energy E, which only contains the 
contributions of the microstructure. The precise formula for E is given in (2.3) below. 
This is done as in |SS03] via an algebraic relation which involves the operator Da dehned 
as follows: 

( 2 . 2 ) 


V^U ■- (V^!l)t+2 - !(V_4'u)t+i , 








where components are understood cyclically (i.e., = 0 ^+ 3 ). In particular, 

V\u = {d2U — iA2u) — i{diu — iAiu) = (V au )2 — au)i ■ 

The operator Va corresponds to a “creation operator” for a magnetic Laplacian in quantum 
mechanics. It was introduced in the context of Ginzburg-Landau by Bogomoln’yi to prove 
the self-duality of the Ginzburg-Landau functional at k = l/\/2, cf. e.g |JT 8 nj . His proof 
relies on identities similar to the next one. 


Lemma 2.1. With the notation above, one has 

~ + P-®3 + X 

and, for any k = 1,2, 3, 

= I^^A«r + pBk + (V X jA)k . 

Proof. We only prove the first relation, the other follows by relabeling coordinates. Notice 
that 


Vau = V-au . 


We compute 


= |(Vam)2 - *(Vam)i| 

= ((Vam)2 - ^(Vam)i) ((Vam)2 + *(Vam)i) 


= |(Vam)i| + |(Vam)2| - *(Vam)i(Vam)2 + *(Vam)2(Vam)i 

= |V>p-z(V»x(^) 


and 


V' X {—iuV'Au) = —iiy'u) X V'aU — iuV' x (V'au) 

= —iiy'u + iA'u) X y'Au) — A'it x y'Au) — hV' x [A'u) 

= -i(V>) X (V» - Wit X (V'm) - \u\y xA'- y'u) X {A'u) 

= -i(V>) X (V» - \u\y X A. 

Since the vector product is antisymmetric, the last expression is real. Recalling the defini¬ 
tion of jA from (1.3) we obtain 


V' X = ReV' X [—iuV'Au\ = —iyy) x (Vam) — |u|^V' x A. 
Adding terms concludes the proof. 


□ 


We now show that the flux of R 3 over every section is constant, due to the divergence-free 
condition. 

Lemma 2.2. Let B G be Qi-periodic and obey divR = 0. Then the quantity 

ip{z) := [ B^dx' 

Jq{z) 














does not depend on z. In particular, if B — bext^s G (implied by the finiteness 

of Eqi^) then for all 2 ; G M, 


/ B 3 dx' = bextL^- 
Jq{z) 

Proof. By the periodicity condition we can test the relation div 5 = 0 (which is true since 
B = V X A) with a function 6 G depending on 0:3 alone. This yields that for any 

9 G Ci(M) 

0 = [ B 39 '{x 3 )dx= f ip{z)9\z) dz . 

It follows that if is constant. The second assertion follows from the fact that B^ dx' > 
p{z)VLfi ' □ 

At this point we are ready to separate the bulk term, and define the microstructure 
functional we shall study below. 


Lemma 2.3. For every admisible pair {u, A) and any parameter set obeying (2.1) one has 

5gl[m, = {nV2bext - + E[u, A], 

where 


E[u,A]\= / {1 — K\/2)\V'jiu\^ + kV2\V\u\^ + \d 3 U — iA^u^ 


' Ql' 


dx 


K 


(2.3) + j \B'\‘^dx + I |53 — fcextpda:, 

jQiXM Jqlx(M.\{0,T)) 

and as above B = V x A, p = \u\‘^. 


Proof. Lemma 2.1 implies 

= (1 — k,\/ 2 )\V'j^u\‘^ + kV^IV^uI"^ + nV^pB^ + k-\/ 2 V' x j'^ . 


The last term integrates to zero by the periodicity of Ja- By Lemma |2.2[ the average of 
the normal component of the magnetic held B 3 is 6 ext- Therefore for each hxed z, we have 


/ IV'Aul^dx' = / 

'< 3 ( 2 ) JQ{^) 


(1 - k.V2)\V'au\‘^ + Ky/2\'D\u\‘^ + K,V2{p - 1)53 + KV^&ext 


dx'. 









We substitute and obtain, using — bext)'^dx = - blJdx, 


Egl[u,A] = 


' Ql,t 


'Ql,' 


(1 - Ky2)|V>|" + kV2|B»m| 2 + tiJj - ^(1 - p)^ 

Idsu - iA3u\‘^ - + fi:\/26e 

|(VxA)fdx+ [ 


dx 


dx 

|(V X A)^ - be^tl'^dx. 


JQlxK JQiX(R\(0,T)) 

Thus, the bulk energy is K\/2bext — the result follows. 


□ 


2.3. Construction of test functions. One main ingredient of the proof of the lower 
bound is the following concentration lemma. This can be seen as a combination of trun¬ 
cation and mollification, and is closely related to Lemma 3.1 from |CCKO08] and Lemma 
2.1 from |CNQ06j . We formulate this lemma in generic dimension, with = {0,L)^, for 
Q2”P6i'iodic functions. In the following only the n = 2 case is used. For / G iloo(l»”) and 
r > 0, we define fr G by averaging over r-balls, 

(2-4) / fiy)(^y- 

I-Dr I JBrix) 

Notice that this operation preserves periodicity. 

Lemma 2.4. Letx G Q'^-periodic, y > 0, 0 < £ < r. Then there is-ilj & 

Q^-periodic, such that 

(i) i){x) > xt{x); 

(ii) sup V’ < sup y; 

r 2nj,n r 

(iii) / ijdy^ —— / y dy; 

JqI Jqi 

(iv) rsup IV-^I < supy; 

(v) r / \\/i/j\dy <-^ xdy. 

Jqi Jqi 

For future reference we remark that these estimates immediately imply 
(2-5) ||^||i 2 (Q«) < j^\\x\\l^{QI)\\x\\l^{QI) and ||VV’||i 2 (Qn) < -^\\x\\L-^iQi)\\x\\LXQl) ■ 


Proof. By homogeneity, it suffices to consider the case supy = 1. Define 


:= |^min|^y2p,l|^ (x) = 


\B. 


r\ JBr(x) 


■ 2 ^-^^ 

min<; ——X 2 riy)A }dy, 


where we use the notation of (2.4). Clearly < 1, and (ii) follows, (iii) is immediate. 

















To prove (i), observe first that xt ^ supx = 1- Fix some x. Since Bi[x) C B 2 r{z) for 
all 2 G Br{x), we have 


Therefore 


Xiix) = 
< 

Xi{x) < min 


BX 

/ xdy 

! Bi{x) 

h. 

/ xdy 

I B2r{z) 


X2r(^),l| 

£n 


£n 


Averaging over Br{x) we obtain 

/ N 1 

Xe{x) < 


\Br 


mm 


' Br{x) 


‘2nj,r. 


\/z G BJx) 


X 2 riz), l\dz = i){x) , 


which proves (i). 

Further, for any pair x, 2 ; one has, writing Cj,(x, z) := {Br{x) \ Br{x + z)) U {Br{x + ;2) \ 
Br{x)), 

h/’ X -%l}[x + z)\< -—-< — Z , 

\Br\ T 

which implies (iv). To prove (v), write analogously 

1 f 


(x) — '0(x + ^)| < 


IB. 


r-l JCr{x,z) 


fdy 


where / = min{2”'r"'f' '^X 2 ri^}- Integrating in x and estimating as above \Cr{x,z)\ 
r'^~^\z\ we get 


< 

r\j 


\'ll){x) — ^^{x + z)\dx < [ < — [ f dx 


'Ql 


>Ql 


Br 


>Q^j 


and the proof is concluded. 


□ 


2.4. The sharp-interface functional. In closing this preliminary section, we introduce 
the sharp-interface version of the functional E and the corresponding function spaces. In 
the sharp-interface functional, a function denoted x represents the characteristic function 
of the normal phase, and is constrained to take values in {0,1}. Thus y is formally the 
equivalent of 1 — p and the approximate Meissner effect (1.1) is imposed via 

B(l-y) = 0. 


Definition 2.5. We say that a pair B G y G BV\oc(^^] {0,1}) is admissible 

for the sharp-interface functional if both of them are Qi-periodic and 

(2.6) divB = 0 distributionally, and B(1 — y) =0 a.e. 

The condition divB = 0 zs understood as B-V6 dx = 0 for all test functions 9 G 












Given an admissible pair iXiB) we set 

(2.7) 




' Ql,t 


Ql,t 




dx + 


'QiX[E\(0,T)] 


\B - be^t^sl^dx. 


This is the sharp-interface functional studied in |CKOn41 ICCKOO^ . In comparing with 
those papers, it is important to notice that several quantities are scaled differently. In 
particular, lengths are rescaled by a factor of L, and magnetic helds by k/\/2. Precisely, 
denoting by x*, B*, E*, e*, L* and 6* the objects used in |CKO04t ICCKO08] . one has 


X*{x) = 1 - x{xL ), 


BGx) = —BixL) 

K 


and correspondingly 


E* = 


2E 


* 2 T 72, 

p* — _ T* — _ h* — — _ h 


3. Lower bound 


3.1. Lower bound with sharp interfaces. In order to prepare some intermediate results 
and to explain the strategy of the proof in a simpler context, we hrst prove the lower 
bound for the sharp-interface functional, recovering the result from |CCKO08] . In proving 
the lower bound the fields y and B will be fixed admissible functions, in the sense of 
Definition 2.5, and we shall simply denote by E the total energy, and by E{z) the part of 
the energy localized in the surface Q{z), i.e.. 


E{z) := / k\D'x\ + 


'Q(z) 


’Qiz) 


Bf + x(B: 


Y 

V2j 


dx'. 


Clearly E{z)dz < E[x,B]. In several lemmas we shall additionally focus on the case 
that the energy is bounded by 


(3.1) 


Flx,B] < -rcb^L^T, 


and that a. z & (0, L) is given, so that 

(3.2) E{z) < • 

The key strategy is to select a good section Q{z) which has small energy. Since the energy 
is small, and the flux of the magnetic held is the same on every section, the magnetic held 
necessarily concentrates on a small subset, whose perimeter is controlled by the energy 
(interior estimate). At the same time, close to the surface the energy favours a uniformly 
distributed magnetic held (exterior term). But “moving around” the magnetic held as ^ 
changes is only possible, due to the divergence-free condition, if the tangential components 
B' are nonzero, which are also penalized by the energy (transport term). Making these 
three ehects quantitative, and balancing them, leads to the lower bound. 


















3.1.1. Equidistribution of the phases. We first show that the average value of x across 
“good” sections is the one that relaxation theory would predict, up to a factor. 


Lemma 3.1. If the admissible pair (x^B) obeys (3.1) and (3.2) for some z G (0,T), then 
(3.3) / 


and 

(3.4) 


' Ql,t 


'Q{z) 


xdx r-u 

be.tL'^T 

K 

xdx' ~ 

bextL^ 

K, 


Proof. We start with the second assertion. Recalling that = X) Jq{z)^Bz — bext)dx' = 0 


(Lemma 2.2) and B{1 — x) = 0, we compute 


'Q(z) 


K 


) dx' 


< 


'Q{z) 


'Q(z) 


’Qiz) 




\ 1/2 
Xdx'j 


'Q{z) 


K 


X ^ - ^3 dx' 


72 


1/2 


K 




72^ 


dl=xdx'j { — I X ( ^ - Rg ) dx' 


K 


'Q{z) 


K 


72 


1/2 


1 

< - 
- 4 


<Q(z) V2 


+ —no 

K 


Therefore, recalling (3.2), 


and 


3 f K , , , ^2 72 , 5, ^2 

7 / H- F{z) < -b^xtL 

4 Jq(,) V 2 k a 


7 / ~^xdx' > bextL"^ - —F{z) > ^bextL 
4 jQiz) 72 K A 


This concludes the proof of (3.4). The proof of (3.3) is analogous, integrating over Ql 
instead of Q{z). 

3.1.2. Interior term. We show that on sections with small energy the magnetic held nec¬ 
essarily concentrates, as captnred by the test fnnction fj. 


Lemma 3.2. For any r > i > 0 and z G (0,T), if holds and is the function 

2\ 1/2 


constructed via Lemma 2.4 from x(-, z) then 


K, 


—= xdx' 
72 jQiz) 

for some universal c > 0. 


'Q{z) 


Bsi^dx'<iF{z) + 


r^ 




K 


FF\z) 


□ R 



















Proof. We write 


K 


Bs = Bsx = —^X + ^83 - 7= j X 


V 2 


K 




Testing with we get 




dx' = —jzz I x'f’ dx' + 
Q{z) 


jQ{z) 

Using Lemma 2.4[i) we obtain 


< 3 ( 2 ) 


K 


B^i- ^]X'^ dx . 


'Q{z) 


xfjdx'P / xXidx' = / x'^dx'p xixe - x)dx' 


'Q(z) 


’QP) 


'Q(z) 


> / xdx' - lx - X^l dx' 


'Q(z) 


'Q(z) 


and therefore 


K 


K 


xdx'< Bsfjdx'p^ / \x-Xi\dx' 


V 2 


'Q{z) 


'Q(z) 


V 2 . 


'Q(z) 


'Q(z) 


K 


^3 - ^ ) xV' dx'. 


The second term can be estimated by 


and the last one by 


'QU) 


'Qiz) 


K 


\X-Xe\dx'<e \D'x\<-F{z), 


'Q{z) 


K 


B 3 - ^]xi^ dx' 


< IIV’IIl 2 (Ql) 


K 




L 2 (Q(^)) 


Collecting terms we obtain 


— I xdx' < I B^ipdx' + IF{z) + \\f)\\L2^^Q^)F^/‘^{z). 


'Q(z) 


'Q{z) 


Using (2.5) and (3.4), we also have 
(3.5) 


iiv’IU 2 (q^) < < (^ 


r 2 


K 


This concludes the proof. 


□ 


3.1.3. Transport term. We now relate the value of B^ over different sections by exploiting 
the f \B'\‘^dx term in the energy. Since we have divi? = 0 we may write 

d^B^ + div'B' = 0 

i.e. B' can be seen as the flux transporting B 3 . Since B 3 takes, approximately, only 
the two values 0 and k/a/2, up to a factor we can understand B' as the velocity with 
which S 3 is transported. Thus we call j \B'\^dx the transport term by analogy with the 
Benamou-Brenier formula for the Wasserstein distance in optimal transport. 














Lemma 3.3. Let B G x (0,T);M^), Qi-periodic, such that divi? = 0. Then for 

any Zi,Z 2 G [0,T] one has 


(3.6) 


{B3{-,Zi) - B^{-,Z2)) fjdx' 


' Ql 


< liv^lh 


\B'\ dx 


'QlX{zi,Z2) 


for any fj G Qi-periodic. The values of B^ on the sections are understood as 

traces, which exist since div i? = 0. 

Proof. This is the same as |CCKOO^ Lemma 2.2], for completeness we give here the short 
argument. We can assume without loss of generality that Zi < Z 2 . We compute, using 
div B = 0 and the Qi-periodicity of B and fj, 


[Bs{z 2 ) - Bsizi)]^; dx' = 


' Ql 


'QiX(^i,22) 


9X2, 


f) dx = 


'QlX{zi,Z2) 


V' ■ B'f; dx 


'QlX(2i,22) 


B' ■V'lpdx < WV'iPWl^^q^) 


\B'\ dx. 


'Qlx(zi,Z2) 


This, together with the same estimate with —f) in place of ip, concludes the proof. □ 

Since B' is nonzero only in a small part of the volume, the embedding of into gives 
an additional factor proportional to 6ext- 


Lemma 3.4. Let (x, 5) be an admissible pair which fulfills (3.1), z G (0,T) such that 
(3.2) holds. Then the function tp from Lemma 3.2 fulfills 


'QP) 


B 2 'ip dx' 


'Q(o) 


B 2 'ip dx' 


< - 


1 




Proof. By Lemma [3.3| we have 
' B 2 'ip dx' — 


We estimate 


’QP) 


\B'\dx = 


'Q(o) 


B 2 'ip dx 


Ql,t 


' Ql.t 


x\B'\dx < 


< 

rsj 


^ ||V'0||l°° 

/ x*) 

•1 Ql,t ) 


\B'\ dx. 


' Ql,t 


1/2 


dx 


' Ql,t 


K 


■) 




To conclude the proof it suffices to recall that ^ 1/r. 


□ 


3.1.4. Exterior term. Finally, we show that the energy outside the sample penalizes con¬ 
figurations with a magnetic field that oscillates strongly at the boundary, as measured by 
the norm. 

















Lemma 3.5. For all admissible and any 'tp G QL-periodic, it holds that 


'Q(o) 


(B 3 - 6ext)V’ dx' 


< 




If (3.2) holds for some z G (0,T) and r > i > 0 and xjj are as in Lemma 3.2, then 


'Q(o) 


{B 3 - be^t)i> dx' 


< 






Proof. Let ip G be Ql periodic, such that V’(‘)0) = 'ip va the sense of traces 

and ||'0||ji^i/2(q^) = IIV'?/i||L2(Q^x(-oo,o)) (we recall that we are using the homogeneous 
norm). Since div B = 0 on (—cx),0) x Ql, using periodicity we obtain 

[ {B - Bext) ■\Ipjdx= [ {B - 5ext)'0 '63 dx' . 

'f<9z,x(—«3i0) d Q{0) 


Therefore 

[ {B 3 - 6 ext)^ dx' 

'Q(o) 

This proves the hrst assertion. It remains to estimate ||t/’||_H'i/ 2 (Q^). By interpolation and 
(2.5) we have 


< ||i? - i?ext||L2(Q^x(-oo,0))||V'?/'||L2(Q^x(-oo,0)) < -^^^^l|V'llr/i/2(Q^) • 




< 


1/2 


, 1/2 


< 


Recalling (3.5) we conclude 


This concludes the proof. 


L2(Q^)||V'0||i2(Q^) 

1/2 


2 

LHQl) 


1/4 


llV^llr^i/2(Q 


< 


Pk 


□ 


3.1.5. Derivation of the lower bound. From what precedes, we deduce the lower bound 
result for the sharp-interface functional. 

Theorem 3.6. For any admissible pair (x,i?) (in the sense of Definition \2. 5| ) and for all 
5ext, iz,L,T > 0 such that 

8 &ext < ^ - o,nd kT > 1, 


F[X,R] > min|6ext«'/'T3AL^&,'/t\2/3Ti/3^2| _ 


one has 


Proof. Let F = F[x,B]. If F > 6ext^7"F^/8 then, since kT > 1, we obtain F > 
bext{f^T)^Bfg and the proof is concluded. We can therefore assume that (3.1) holds. 
By a mean-value argument, we may choose z G (0,T) such that F{z) < F/T, so that 
(3.2[) holds as well. We start by constructing xp as in Lemma 3.2, namely, with Lemma 2.4 






















applied to some 0 < £ < r chosen below. One key observation is that, since "0 

depends only on x', 


'Q{z) 


B'i'lp dx' = / &ext'0 dx' + 


' Ql 


/ B^'ip dx — / B^ip dx' 
Jq{z) JQiO) 

With Lemma 13.41 and Lemma 13.51 this can be tranformed into 

1 


+ / {B 3 - bext)^ dx'. 

Jq{ 0 ) 


^ B:d’ dx' < jl dx' + ci + c F ^!^, 


Lemma 2.4[iii) gives 


' Ql 


bext'if^dx' < bextMiXQL) < 4—6ext / X dx'. 

t Jq{z) 


Recalling Lemma 3.2 


/ Xdx' < [ Bsijdx'+ c£F{z)+ F^/^{z). 

V2jQiz) Jq(z) V / 


Combining the last three estimates gives 


K 

71 


'Q{z) 


xdx' <4-—bext / xdx' + c- 


1 


£2 


'Q{Q 


K 


FF^ + 




£2 


K, 


i?l/2 


+ c£F{z) + c 


r2 6..W2xi/2 


£2 K 


F^/^{z). 


Assume now that 
(3.7) 


0 < £ < r < 


K 


8 h 


ext 


1/2 


so that the coefficient of X dx' in the right-hand side is smaller than the one on the 
left-hand side. This is possible, since we assumed 6ext < ^/8. Then, recalling (3.4) and 
F{z) < F/T, 


b^xtL^ <- 


1 / 6 ..W 2 yx 1/2 




^1/2 ^ 


1/2 , .F , /r2 6exti:'\'/2FV2 


£2 


K 


fF^ + £- + 

' rri ' 


£2 


K 


yl/2 • 


At least one of the terms in the right-hand side has to be at least one-quarter of the one 

ii( 

£^k 


on the left, and therefore for all pairs {r,£) which obey (3.7), we have 

,2 




F > min bextL —, bextL —, b^^tL —, bextL 


F > nbextL'^T min 


r2 £2 1 £2>. 


Equivalently, 

(3.8) 
























We finally have to choose r and i, and check that in each case some terms give the optimal 
bound, and the others are irrelevant. Balancing the hrst three terms we obtain 


This choice is admissible only if (3.7) is satisfied, which since kT > 1 is equivalent to 
< {K/SbextY^'^- In this case, (3.8) becomes 

F > K^/\xtL‘^T^^''mm{l, 1,1, (kT)^/^} 

and since kT > 1 the assertion holds. 


If instead > (n/SbextY^'^, we choose r = i{K/8bextY^‘^- Inserting this into (3.8) 

and then balancing the first and third term yields, after some rearrangement. 


F > min 11, {kTY'^ ' , 1, {kT) 


2/3 


bext 

K, 


4/3' 


We observe that since kT > 1 and (kTYF > {n/SbextY^"^ terms of the form 


(/s:T)“ ( — 

K 


with a > “1(3 are bounded from below. This concludes the proof. 


□ 


3.2. Lower bound for the Ginzburg-Landau functional. The proof is structured in 
a similar way as the one for the sharp-interface functional, but contains several additional 
difficulties. In particular, when working with diffuse interfaces we can enforce the Meissner 
condition only in a weak sense, see Section |3.2.2[ This generates difficulties both in the 


interior estimate and in the transport term. Additionally, the energy does not directly 
control the size of the boundary of the normal phase, and a suitable estimate needs to be 


formulated and proven (Section 3.2.3). The only term which can be treated in the same 
way is the one corresponding to the external field. 

We consider a pair (u. A) G with Eqi\u,A\ finite, hence E[u,A\ finite. By Lemma 
3.7 below, we can assume that p < 1 pointwise. For any G (0, T), we denote by E{z) the 


energy contained in the section Q{z) := Ql x {z}, 


E{z) := 


’QE) 


(1 - kV2)\Vau\‘^ + kV2\V\u\^ + ( B: 




p) 


dx' 


' Q{z) 


\\B'\^ + \d-^u — dx'. 


We write for brevity E = E[u,A\. We recall that p = \u\^ and dehne 
(3.9) x--=il-pf- 















3.2.1. Normalization of the density. We first show that we can assume without loss of 
generality that |m| < 1 , or, equivalently, p < 1 . 


Lemma 3.7. Let {u,A) be an admissible pair. Then, 



|V(7''")l<|V^ti|, 

VI 

w 

Vam 

Let u : 

X (0,T) C be defined by 




1 u{x) 

if \u\{x 

)<1 

(3.10) 

“ (;r) 

U«l 

if \u\{x 

)> 1 . 


Then {u,A) is in and E[u,A] < E[u,A]. 


Proof. By density it suffices to prove all assertions under the additional assumption that 
M G at points where p 7 ^ 0. Writing locally u = one obtains 

Vau = [Vp^/2 + ip^/2(V0 - A)] , 

hence locally 


(3.11) \VAu\^ = \Vp^/^\^ + p\Ve-A\\ 

The second assertion follows immediately from the dehnition of Va in (2.2). It also follows 
directly from (3.11) that \Vau\ < |VA^i|, and hence Egi[u,A] < Egi^[u,A]. Since E and 
£^gl differ by a constant, we conclude that E[u, A] < E[u, A]. □ 


3.2.2. Meissner effect “on average”. In the reduced model we had the compatibility con¬ 
dition B{1 — x) = 0. In the true model, we expect as in (1.1) that Bp is, in some sense, 
small. The next lemmas make this quantitative, in appropriate weak norms. 


Lemma 3.8. Let {u,A) be admissible, and assume (2.1) and p < 1. For any z, we have 


(3.12) 

and 

(3.13) 


'Q{z) 


pB^ dx' 


< SE{z) 


'QiP 


pB^il - p)dx' 


< 16 E{z). 


Further, for any p G L°° fl x (0,T)), Qi-periodic and such that p = 0 on Q{0) 

and Q(T) (in the sense of traces), and any k = 1,2,3, we have 


(3.14) 


pBkp dx 


' Ql,' 


< 




Proof. We start with (3.12). By Lemma [2T| 


'Q(z) 


pBs dx' = / [I V>r - V' X - \V\u\^] dx'. 


’QF) 















The integral of V' xJ'a is zero, since is Qi-periodic and only in-plane derivatives appear. 
The other terms are bonnded by the energy, i.e.. 


1 

2 




\V\u\‘^dx' < 


Iq(z) - 1 -kV 2 


< mz) , 


where we nsed k < 1/2. This conclndes the proof of (3.12). 
The argnment for (|3.13|) is similar. We write 


'Qiz) 


pB,{l - p) dx' = / - V' X J'a - (1 - P) dx' . 


'Q{z) 


The terms with Vau and TPau can be estimated as above. The term with V' x however 
needs more care. Since Ja and p are Q^^-periodic, IjaI < au\-, P < 1 and < 

|Vam|, an integration by parts leads to 


'Q{z) 


(1 - p)V X j'^dx 


'Q(z) 


X dx' 


< 


K 


x/2 


E{z) 


This proves (3.13). 

Finally, we consider (3.14). Here we inclnde the other components, and we reqnired that 
the localization fnnction p vanishes on the top and bottom bonndaries (where we have no 
periodicity). We compnte 


/ pBkp dx 

— 

/ 

- \vfu\'^ - (V x jA)k 

p dx 

J Ql,t 


■IQl.T 




< 16 E||(p||l°o + f IjAllVpIdx 

J Ql,t 

< 16 E\\p\\loo + ■ 


□ 


3.2.3. Surface energy. We now show how the snrface energy can be recovered from the 
fnnctional. This arises from the combination of a and a (1 — p)^ term, bnt the 

latter is not directly present in the energy, and needs hrst to be reconstrncted. 


Lemma 3.9. Let (n, H) he admissible, and assume (2.1) and p < 1. For every z G (0,T) 
we have 


(3.15) 


/ p{l-pfdx'<E{z). 
Jq{z) 


The function x defined in (3.9) satisfies 


(3.16) 



Vx\dx'<E{z) 






















and, for any i > 0, 


(3.17) 


/ Ix^ - Xl dx < -E{z). 

iQiz) 


Here xe is the average of x over i-balls, defined in (2.4). 
Proof. We have 






^3-^(1-p) > P\B3 -^(1-p) 




72' 


> yP(l - P)^ - kV2pBs{1 - p). 


Therefore, recalling (3.13), we have 


K 


— / p{l — pY dx < E[z) + 1Q\/2 kE[z) , 


’QP) 


and (3.15) is proven. 

At the same time, since < |Vam|, we have 


’QP) 


\Vp^/^\‘^dP <4:E{z). 


Therefore 


E(z) > / [\Vp'/X + K^p(l-pf]dx' 


'Q{z) 


> 2 / Kp^^‘^{l — p)\'V'p^^‘^\dx' = I iV'xldx', 


'Q(z) 


'Q{z) 


since V'x = — l)V'p^/^. This proves (3.16). 

Let now xe be defined as in (2.4). By Jensen’s inequality and the mean-value theorem 
we obtain 


'Q{z) 


\Xe - x\ dx' < sup / \x{x' + h') - x{x')\dx' <i |V'x| dx'. 


\h'\<eJQ{z) 


'Q(z) 


Notice that this still holds if 7 > L (indeed, in this case the coefficient could be improved 
to L). This concludes the proof of (3.17). □ 


3.2.4. Equidistrihution of the phases. We show that in every section Q{z) with a good 
energy bound the volume fraction of the normal phase is approximately “right”, in the 
sense that it can be obtained assuming that Us equals k/\/2 in the normal phase, and zero 
outside. 


Lemma 3.10. Let {u,A) he admissible, and assume (2.1) and p < 1. 









(i) If 


(3.18) 


E{z) < -Khey±L^ , 

o 


then 


xdx'= I {1 — pY dx'^L‘^ 


>Q(z) JQ{z) 


K 


(ii) If 


(3.19) 


E < -nbextL'^T , 


then 


xdx = 


Ql,t 


Ql,t 


K 


Proof. We start with assertion (i). Using Lemma 2.2 we write 


h ^ 

OextJ^ -^ 


72 


[ xdx'=l (B,-^{l-pAdx' 

'q(z) Jq{z) V V2 / 


’QP) 


K 


(1 - p) UBg-^(1 - p) dx' + / pEg dx' 


72 


'QP) 


and estimate, using Cauchy-Schwarz, (3.12), and k < 1/2, 

\ 1/2 


- A 


72 


'QP) 


xdx' 


< 


([ (l-p)2da;') E^/\z) + 8E{z) 

\Jq{z) j 




Q{z) 


^ , E{z) E{z) 
xdx H-^ +4^^ 

K K 


The conclusion follows. 

The second part is proven analogously, just extending all integrals to Ql,t- 


□ 


3.2.5. Interior term. The next lemma is the key ingredient of our proof. It shows that, if 
the localization is performed appropriately, in low-energy sections the held necessarily 
concentrates. The concentration is made quantitative by a test function constructed via 
Lemma 2.4 starting from y = (1 — p)^. 


Lemma 3.11. For any admissible {u, A) with p < 1, any parameters which obey (2.1), 
any r > £ > 0, and any z G (0,T) such that (3.18) holds one has 

72 Jq{z) Jq(z) ^ \ I J 

The function fj is the one is obtained via Lemma 2.4 from the restriction to Q{z) of 
x = {i-P?- 














Proof. We first write 


K xdx' = K (x - X^)dx' + K xix- Xi)dx' + k xxdx '. 
jQ{z) Jq{z) Jq(z) Jq{z) 

We now estimate the three terms on the right-hand side. For the hrst one, we compute 
[ ix-X^)dx' < k[ {1-pf{2p-p‘^)dx' 

jQ{z) JQ{z) 

^ 2 < dd{z) 


<2k p{l — pY dx‘ 

Jq{z) 


K 


by (3.15). For the second, 


xix- xY)dx <K lx - Xe\dx' < iE{z) 
jQ(z) Jq{z) 


by (3.17). For the third, we write, recalling Lemma[2.4Ki) and the dehnition of x and Xii 


Therefore 
(3.20) 

At this point we write 


^ / XXedx' <K (1 — p)'il)dx' . 
Jq{z) Jq{z) 


K f X dx' < K f (1 — p)'ijj dx' + c—^ -|- ciE{z). 
Jq{z) Jq(z) 


[ Bs'iljdx'=f -^{1 - p)Y dx'+ [ (B 3 --^{1 - p)\Y dx'. 

'Q(z) JQiz) v2 Jq{z) \ X2 J 


We use 


'Q(d 


(^3 - ^(1 - Ydx' < \\Y\\LXQr.)E^^^iz 


and (3.20) to obtain 


/ Xdx'< [ BsYdx'+ + ciE{z)+ \\Y\\mQL)E^^^iz). 

X2Jq{z) Jq{z) 




Recalling ( |2.5[ ), 

J.2 

ll^lli^Qi) ~ ~p\\x\\lYQl) < ^ 2 ^ ’ 

where we used once again Lemma 3.10Ki), concludes the proof. 


□ 














3.2.6. Transport term. It remains to relate the behavior of in the interior with the 
behavior at the boundary. As in the sharp-interface case, this is done in two steps, but 
both steps are different than the corresponding ones in Section 3.1.3 The hrst estimate is 
easy, but does not give the optimal bound, since it is oblivious to the fact that B' needs 
to be concentrated on a small volume (this is relevant, since we are estimating an term 


with an energy that contains the corresponding 
in the following lemma. 


norm 


The estimate is then improved 


Lemma 3.12. Let {u,A) be admissible, -0 G Qi-periodic. For any pair zi,Z 2 G 

(0, T) we have 


'Q(2i) 


B^tl: dx — 


'< 9 ( 22 ) 


B^%Ij dx' 




Proof. We can assume without loss of generality that Zi < Z 2 . We compute, using div B = 0 
and the Q^-periodicity of B and ip, 


' Ql 


[BsP, Z 2 ) - BsP, zi)]f) dx' = 


[ ^ 

Qlx{zi,Z2) ^^3 


Ip dx = — 


' QlX{z-i_,Z2) 


V ■ B'fj dx 


'QlX(zi,Z2) 


B' ■ Vip dx 


< 


'QlX{zi,Z2) 


\l/2 / p \ 1/2 

B'l^dxj 1^12^2 — 2 :i| J I'Vipl'^dx'j 


This concludes the proof. 


□ 


Lemma 3.13. Let {u, A) be admissible, and assume (2.1) holds and p < 1. If (3.19) holds, 
then for any z G (0,T) and any tp G Qi-periodic, we have 


'Qiz) 


B^pj dx' 


'Q(o) 


B^^p dx' 


< 




K, 


\V^P\\L^E^/^ + \\V^P\\L^Q,)E^/\ 


Proof. If z < 2 , this follows immediately from Lemma 3.12 Assume z > 2 , and hx 
6 G (0, z/2). Let 77 : M —>■ M be dehned by 


Vixs) := < 


5 

1 

z - 


X 3 


if 0 < xs < (5, 
if 6 < X 3 < z — S, 
if z — 6 < X 3 < z, 
otherwise. 


0 















We compute 


' Ql,t 


B^{x)'iIj{x')^{x^) dx = - 


[ 1^' 

'Ql,t 9x3 


x)'ilj{x')rj{x^) dx 


/ V' ■ dx = — B' ■ (VV)^ dx 

'Ql.T JQl,T 

-f pB'■ {Vij)r]dx - f {1-p)B'■ {V^p)pdx. 

•^Ql,T '^Ql,T 


The second term can be estimated by 


Ql,t 


(1 — p)B' ■ {y''ip)pdx 


r 

<1 / (1 — p)^ da: I 

'Ql,t j \JQl,t 


Since (3.19) holds, Lemma 3.10 gives 


/ {l-pfdx< 

'Ql,t ^ 


\ 1/2 

B'l^dx] IIVVIU- 


K..tL^T 


and therefore 


(1 — p)B' ■ {V''4))pdx 


' Ql' 


h r 271 \ 1/2 

< ( -j 


K 


For the hrst term we use Lemma 12.11 to obtain 


pB’ ■ {V''ip)pdx 


' Ql: 


< / {\yAu\^ + \VAu\^)\V'lP\pdx 

d Ql,t 


(V X Ja) ■ {V'ip)pdx 


' Ql: 


The part containing ja can be transformed according to 


/ p (V X Ja) ■ VV dx = — Vp ■ Ja x VV dx. 

JQl.t JQl,t 

This is easily proven by integration by parts and using the symmetry of the triple product. 
Since ||V> = \\V'4 )\\l2 (q^)\\^^\\l^P3,t)) and ||jA||L2(Qi,T) ^ we obtain 


pB' ■ {y''tp)p dx 


'Ql: 


< 


£;||vV>lk- + £‘'"l|vVllL.,Q.)ll^lli»((o.T)). 



















Adding terms, and replacing rj by — 77 , we conclude that 

'h„,eTV'^ 


(3.21) 


1 dx?. 


< 

r\j 


] £''^||VVIU« + b''''I|vV>IU.(o.)II^IIi.((o,t)) 


K 


where we dropped the term i?||VVII l^ using k < 1 and the assumption (3.19) and we 
dehned 


fixs) := 


’Qi^s) 


B'i'ip dx' 


This estimate controls the variation of / on a scale 5, and can be combined with Lemma 


3.12, which gives a bound on the Holder 1/2 norm of /, to obtain a pointwise estimate. 


Precisely, since 

we get 

Therefore 


f-p-dx3 = \ [ fdx3-\ [ fdX3, 


dx 


3 dx 3 


z-5 


+ sup |/(x 3 ) -/( 0 )| + sup |/(a; 3 ) -/(^)|. 

3:36(0,( 5 ) 2:36(2—5,2) 


i/(=) - /(o)i < ^ £‘'"iivviu» + 


We hnally choose 5 = 1 and conclude the proof. 


□ 


3.2.7. Exterior term. 

Lemma 3.14. For all admissible {u,A) we have 



be^t)' 4 ^dx' 




Proof. This is the same as in the sharp-interface case, cf. Lemma 3.5 


3.2.8. Proof of the lower bound. 

Theorem 3.15. Let (n, A) be an admissible pair. If 

5ext < and kT >1 , 

o Z 


□ 


one has 


E[a,A] > min blit 














Proof. By Lemma 3.7 we can assume without loss of generality that p < 1. We can assume 
B — Bcxt € L‘^{QlX otherwise the energy is inhnite, so that we can use Lemma 


2.2 


If (3.19) does not hold then E > bex_tL‘^{K.TY^'^, since kT > 1, and the proof is concluded. 


Therefore we can assume that (3.19) holds. We choose 2 ; G (0,T) such that E{z) < 


E\u , A]/T, so that in particular (3.18) holds. Let '0 he the function constructed via Lemma 
2.4 from the restriction to Q{z) of x = (1 ~ pY, as in Lemma 3.11, for some parameters 


r, i still to be chosen. We start from the identity 


/ dx' = / fcextV' dx' + 

Jq{z) Jql 

From Lemma [ 3 . 131 


/ dx' — / i? 3'0 dx' 

Iq{z) Jq{o) 


+ {B 3 - 6ext)V' dx' 

Jq{o) 


/ B^ip dx' — / B^f) dx' 
Iq{z) Jq{0 ) 




K J 


+ 


be^ 


2 \ 1/2 


E^/\ 


( 2\1/2 

where we used || VV||l°° ^ l/'i’ and || VV||l 2 (( 3 ^) < ( j (both obtained from Lemma 
|2.4[iv) and (v) and Lemma 3.10|(i)). Analogously, Lemma 3.14 shows that 


'Q(o) 






kP 


where we inserted the bound on the H^P norm of ip : 

imny^Q.) i ii’AiiFq.) iiw-ii7k, :s (■ 

Further, by Lemma 2.4[iii), 


r6^2^i/2 

kP 


' Ql 


I I 

be^tfjdx' < feextllV'llLpQi) < 4&ext^ / X dx' , 


'Q{z) 


therefore 


'Q{z) 


K 1,1 ' <A,. ^ (K.XX'- 

Bslpdx <4&ext^ / xdx +c- 


E^P 


'Q(z) I’ V ^ 

' r r2 \ 1/2 / 7 r2 \ 1/2 

+ C| ^1/2 + ^ / ^^extL \ 


kP 


kP 


E^P. 


At this point we use Lemma 3.11, which states that 


A / Xdx'<! B# dx’ + + c (lAA/ 

v2 Jq{z) Jq{z) ^ \ ^ J 




































Combining the previous estimate gives 


[ Xdx' <Ah^J— [ xdx' 

V2 Jq(z) Jq{z) 


'Q{z) 

Assume now that 
(3.22) 


+ c 


Q(z) 
rbe^tL 

kP‘ 


+ c- 
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2\ 1/2 


K 

E{z 
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2\ 1/2 


EI/2 ^ ^ ^ ^ 


0 < £ < r < 
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be^tL" 

kP 




K& 


^1/2 


86 , 


ext 


1/2 


SO that the hrst term on the right is no larger than the one on the left divided by a/2. 
Recalling Lemma 3.10[ i) and E{z) < E/T we conclude that for all pairs {r,i) which obey 
(3.22), we have 


be.tL^ <- 


1 fb ,,,L^T y^ 


ee^ + 


6ext^^ 


2\ 1/2 


kP 


^1/2 




-^ext 

kP 
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„E 


E^!'^ + - + £- + 
kT T 




kP 


yl/2 • 


We remark that only the second and the fourth term are new with respect to the sharp- 
interface case. At least one of the six terms has to be at least one-sixth of the total, 
therefore 

, r2 rr. . r2^, . 2 1 

-,6ext^ KT,6ext^ -,&ext^ 

r 


E > min <j 6ext^^^, 6ext^^K^^ 6ext^^ —, be^tL'^nT, be^tL^^, 6ext^^’ 


/■ 


Equivalently, 


P P P 1 P 

E > K6extf^"T min { —, -, —, 1, -, - 

The fourth term can be dropped, since the sixth one is always less than 1. Therefore we 
can focus on 

„2 o 2 02 £2 

T2’ T’ rT’ 'PPP 


(3.23) 


E > nbey.tL'^T min <j —, —, -, -P-, /- 


where a dot marks the term we already know to be irrelevant. In comparing with (3.8), 
we see that the only new term is the second one, P/T. Averaging the hrst and the sixth 
we see that ^|T would be irrelevant; hence the second term is irrelevant for all choices of 
£ > 1. We shall see later that this is the case. 

We hnally have to choose r and £, and check that in each case some terms give the 
optimal bound, and the others are irrelevant. Since we already know the scalings, we 
do not need to check all possible combinations. Balancing the hrst, third and hfth term 
suggests the choice 


£ = T^/7«:-3/7^ 


r = T^n^-Pp 



















This choice is admissible only if (3.22) is satished. The hrst condition is always true, since 


kT > 1; the second one is equivalent to Since k < 1 < kT, one can 

compute i = /K > 1, hence the second term can indeed be dropped. In this case, 


(3.23) becomes 


E > «=*/"6exti:'T3/7min {l, •, 1, •, 1, . 

Since kT > 1, in the regime < (/t/Sfeext)^'^^ we have shown E > 


If instead {kTY^ > (^/gfeext)^'^^, we need to choose r = i{K/8bextY^‘^- Then (3.22) is 
always satisfied. Balancing the first and fifth term in (3.23) with this constraint results, 
after some rearrangement, into 




2/3 


and 




mm 
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K K \ K 


7/6 /, \4/3' 

\ K 


Again, a dot marks terms we already know to be irrelevant. We observe that since kT > 1 
and {kTY/'^ > (^/Sfeext)^'^^ one has 


(«T)“ > 1 


whenever a > ‘^(3. This permits to show that the second, third and sixth terms do not 
contribute, and therefore concludes the proof. □ 

4. Upper bound 

Before presenting our construction for the energy Uql we construct fields with optimal 
scaling for the sharp-interface functional F. We prove a refined version of the results of 
[CKOOdj . giving a construction which satisfies several additional properties, which will be 
needed in the following generalization to Uql- In particular, we need to make sure that 
there is an integer number of flux quanta in each flux tube and that the thickening of the 
tubes on the scale of the correlation length still has small volume. The construction for 
Uql will then be derived from this one. 

4.1. Construction with sharp interfaces. The key point in the construction is to use 
in two stages a subdivision of the domain into subsets with integer flux. This is done via 
the domain subdivision algorithm of Lemma 4.2 We first subdivide the domain down to 
scale L/N, and obtain rectangles (rj)j=iwhich will be used to £x the microstructure 
in the central section of the sample, QlX {T/2}. In each of these rectangles, the magnetic 
held will be concentrated in a smaller concentric rectangle fj, keeping the same flux. Then 
we subdivide a second time, again using Lemma |4.2|, down to smaller rectangles, which 


will be the ones used close to the surfaces of the sample, Ql x {0} and Ql x {T}. At the 
same time, the total flux inside each rectangle is concentrated in a smaller rectangle, so 
that the intensity of the magnetic field is the one preferred by the energy, k/\/2 (see Figure 
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Figure 2. Sketch of the notation used in the proof of Theorem 4.1 


[^. After this setup we will make the actual branching construction, which corresponds to 


the subdivision generated in the second application of Lemma 4.2 


Before stating the main result of this section we recall the dehnition of F in (2.7) and 


introduce the notation {u)p for a p-neighbourhood of a set u. Precisely, for u C and 

5 > 0 , 


(4.1) 


(a;)^ := {x G : dist(x, a;) < ^} = Bs{x) 


xGcu 


where the distance is interpreted QL-periodically, dist(x. A) := m.i{\x — z — kL\ : z G A, k E 
X {0}}. 


Theorem 4.1. For any 6ext, k, L,T > 0 such that 


and 

(4.2) 


2^ext < K < L 


L > min 


kT > 1, 


e 27rZ, 


8^2/3 I 


there is a pair {x,B), admissible in the sense of Definition 2.5, and such that 

F[X,5] <min|6ext«'/'T3AF2,6,'/t\2/32^i/3^2| _ 


The set u ■.= {x E Ql,t '■ x{^) = 1} is formed by the union of finitely many sheared 
parallelepipeds, with two faces normal to e^, and 


(4.3) 


(w)i/„\w|<-/ |C\i 

^ JQl,t 


For any z E (0,T) the flux of B across each connected section of {x^ = z}r]u is an integer 
multiple of 271. 


Before giving the proof of Theorem 4T, we formulate and prove the partial results that 
will be needed. We start with the domain subdivision. Here we rehne a flux pattern 






















making sure that each component keeps the quantization condition. In order for the field 
to maintain the optimal intensity, the areas are changed. Some parts of the construction 
would be simpler if one would work with squares, but then the macroscopic distribution of 
the flux would be modified. We work with rectangles, of aspect ratio uniformly close to 1, 
so that on any scale the perturbation to the distribution of flux is kept to a minimum. 


Lemma 4.2. Let i?o,i = (0, a) x (0,6) and > 0 be such that |a < 6 < 3a and obB^ G 
27rZ. Then for any k E N there are 4^ pairwise disjoint rectangles {Rk,i}i=i,..., 4 >=, Rk,i = 
Xk,i + (0,afc,j) X such that \ak^i < bk,i < Sok^i, each R^^i is the union of four Rk+i,i 

(up to null sets), and ak,ibk,iB^ G 27rZ for all k,i. Further, \ak,ihk,i — 4“*^a6| < 4ti/B^. 


Proof. The lemma follows immediately from Lemma 4.3[ if one ignores half of the stages. 
Precisely, if are the rectangles produced by Lemma 

for i = 1,..., 4^ = 2^^. 


4.3 


we set Rk,i = R* 2 k,v 
□ 


For the proof it is more convenient to focus on the following version, in which only one 
subdivision of each rectangle into two is performed at each step. To simplify the notation, 
we say that a rectangle R = x' + (0, a) x (0, h) C is i?*-good, for some B^ > 0, if 

(4.4) -ci<b <3a and abB^, G 27rZ. 

o 

Lemma 4.3. Let B^ > 0, and let i?o,i be a B^-good rectangle. Then for any k E 'N there 
are 2^ pairwise disjoint B^-good rectangles such that each Rk^i is the union 

of two Rk+i,i (up to null sets) and \ \Rk,i\ — 2“^|i?o,i|| < dvr/i?*. 


Proof. The construction is iterative, starting with Consider one rectangle at step 

k, say, Rk,i = Xk,i + {0,ak,i) x (0,6fc,i). Assume for definiteness that ak,i < bk,i. If Rk,i 
is empty, i.e., has side lengths zero, we replace Rk,i by two empty rectangles, setting 
Rk+i, 2 i ■= Rk+i, 2 i+i ■= 0 - If B^ak,ibk,i = 27r, we replace it by a copy of itself and an empty 
rectangle, setting Rk+i^ 2 i '■= Rk,i and Rk+i^ 2 i+i ■= 0 - Otherwise, we set 

Bd/iOk^ibk^i 

dvr 

and replace Rk^i by the two rectangles Xk,i + ( 0 , 0 *,^*) x (0,?/) and Xk,i + ( 0 , 0 ^,*) x {y,bk,i). 
Here [fj = max{z E Z : z <t}. 

Clearly y < bk^i/^. At the same time, since for any x G N with z >2 one has [^/2J > x/3, 
one has y > bk^i/3. Therefore the aspect ratio of the two new rectangles is also not larger 
than 3 and they are i?*-good. 

It remains to estimate the area. Consider one rectangle Rk,i at stage k. It has been 
generated from i?o,i by k subdivision steps. Let Sh be the area of the rectangle at stage 
h along this subdivision path, so that sq := |.Ro,i|- By the dehnition of y, we obtain 
|2sh — s/i-i| < dvr/H*. Summing the series we obtain 


27r 

Bd^ak i 


|2"sfc - |i?o,i|| < |2^/. - 2^-hh-i\ < < 2"=^^ 


h=l 


h=l 


Bd, 












which gives \sk — 2 ^|-Ro,i 


□ 


We next estimate the energy of a flux conhguration on the boundary. We assume that 
the magnetic held inside each rectangle Rj is concentrated in a subrectangle Vj. The 
difference of the helds ajlr^ — then has average zero over the larger rectangle Rj. 

Here and below we denote by 1^; the characteristic function of a set E, 1e{x) = 1 if x E E, 
0 otherwise. 


Lemma 4.4. For M G N let 


{Ri,...,Rm} and {ri,...,rM} 


be reetangles such that each of them has aspect ratio not larger than 3, rj C Rj C Ql, with 
the Rj pairwise disjoint. Let aj,Aj eM. be such that aj\rj\ = Aj\Rj\. Then 


M 


i=i 


2 


M 








Proof. Fix one index j, and let g := Ojlr^ — Ajl^^. Since g has average 0 over Rj, for any 
if E L^{Rj) we have 


/ 9T dx' 


1 

-G 

o 


1 

-G 

o 

— I®il \\T 9^o||Ll(rj-) 

Jr, 


\JRj 


\Jrj 



where (po is the average of p over Rj. Since the trace of a function (in three dimensions) 
belongs to if p E E[l^^{Rj x (0, cx))) we obtain 

\\T - <^o||Li(r,) < Vjf^Wp - < \rjf^Vp\\L2(R^>po,oo)) ■ 

Since the aspect ratio of the rectangles is controlled, the constant is universal. 

Now £x $ G X (0, cxd)), Qi-periodic. Let pj be the average of $ over Rj. Then 

the same computation gives 


JQl ,• Jr. 


aj(<h — pj)dx' 




kil kjf ^1 |v$|U2(k^.x(o,oo)) 


z,x(0,oo)) ) 


where in the last step we used Cauchy-Schwarz. 


□ 


Controlling the norm of the normal component of B on the boundary is sufficient 

to estimate the energetic cost of the magnetic held outside the sample. We recall this 
general fact in the following Lemma. 



Lemma 4.5. Let g G Q^-periodic, with average 6 ext £ Then there is B E 

also Qi-periodic, such that divB = 0, B^{x',x^) = g{x') for x^ > 0, and 


' Qlx(— oo,0) 


\B - be^t^sl^dx 


< 


\\g - &ext 


||2 


Proof. It suffices to consider the case 6ext = 0. Let g{k') be the Fourier coefficients of g, so 
that g{x') = T^k'^o I^P(^')/I^'I ~ \\9\\h-i/2- Since g has average 

0, ^(0) = 0. We dehne, for xs < 0, 


B(x',X3):= P'^'-^'B{k',Xs) 


fc'e27rZ2/L 


where 


k' 


B,{k\ xs) := g{k')e \^'\^^, B\k\ X 3 ) := i^^g{k')e\'^'\^^. 
It is then straightforward to check that the stated properties are satished. 


□ 


Before starting the construction in the interior region, we introduce a separate notation 
for the interior contribution to the energy. We dehne, for fl C open, y G BV{Q) and 
B G L2(fI;M3), 

(4.5) F“\x,B,{l]:^ f k\Dx\+ f ' " 


K 




dx. 


so that F[x, B] = B, Ql^t] + /( 


(9i,x[K\(0,T)] 


\B — 6 exte 3 pdx. By Lemma 


4.5 


it suffices 


to control F”*[x, 5, + \\B - &ext||^-i/ 2 (Q^x{o}) + ~ ^ext||^-i/ 2 (Q^x{r})- 

The next construction step is a procedure to generate an admissible y and B with given 
boundary data in a slab x (0, t). The explicit construction is done for the case that the 
boundary data are characteristic functions of rectangles. An extension to the case where 
circles are used, which gives a smaller surface energy (by a factor which does not affect the 
scaling), is discussed in |CGOS] . 

Lemma 4.6. Let r := p+ (0, a) x (0, b), r := p+ (0, d) x (0, b) C be two rectangles with 
|r| = |f|, a ~ 6 and d ~ 6 . For any t > a there are (y, i?) : x [0,t] —)■ {0,1} x such 

that 


K 


B3(-,-,o) = ^w, 


K 


Bz{-, •, t) — , 


divS = 0 , 

the first two in the sense of traces, and such that, defining uj := {x : B{x) 7 ^ 0}, 




\p — pp + ab 


If R = (xi, yi) X (x 2 , 2 / 2 ) is such that r Uf G R, then u G R x [0, f]. For any S < a. 


with {u:)s defined as in {4-1)- 









Proof. We define 99 : M —)■ M by 


ip^xs) := exp 




Then (^9(0) = 1, (p{t) = d/a = b/b, \ip'\{x^) < 1/t for x^ G [0,f]. Fnrther, we define 

^(Vi-Vi\ ( xMx 3)\ 
v{x) := — I P 2 -P 2 1 + \ X 2 /^{X 3 )\ , 

which is a diffeomorphism of x [0, t] into itself, with det Vn = 1 pointwise, and hnally 
dehne y and B by 

(4.6) x{v{x)) := W(a;'), B^Hx)) := , B'{v{x)) := ^l,{x')d 3 v'{x). 

Let 9 G X (0,t)). By a change of variables 


/m2x(0,P 


{d3eB3 + ve-B')dx = 


/r2x(0,P 


{836 o V B 30 V P V'9 o V ■ B' o v)dx . 


Inserting the dehnition from (4.6) this becomes 


K 


—= / [836 O v\r + '^'9 o V ■ {lrd2,v')]dx = 

\/2 jR2x(0,t) 


/ lr{x')-^{9 o v){x)dx = . 

dX3^ 


Therefore the divergence condition is satisfied. 

Finally, ||i?'||ioo < fi:||lrcl 3 n'||Loo < K,\p — p\/t + cna/t. Since a ~ 6 and the volnme of its 
snpport is tab, we obtain 

f \Bfdx<P’^xil±^ab. 

iR2x(0,p ^ 

From the dehnition of u one easily obtains the other properties. □ 


At this point we present the branching constrnction, which gives the rehnement of the 
magnetic hnx close to the bonndary. We start from a prescribed interior structure, which 
can be obtained either by uniform subdivision of the domain and quantization of the held, 
or by nonuniform subdivision of the domain with a uniform held. We stress that the 
rectangles in which the construction is localized do not need to cover Q^. Indeed, for 
very small helds the optimal scaling is only obtained if these rectangles cover a very small 
fraction of Ql (the volume fraction will be 7 ^ in the proof of Theorem 4.1 below), see 
Figure 


Lemma 4.7. Let k, L,T,do, po, N > 0, with N'^ G N and n <1, and assume that two sets 
of rectangles are given, 


{ri,... ,rAr 2 } and {fi,...,fAr 2 } , 





Ql X {T} 


Ql X {0} 




Figure 3. Sketch of the horizontal geometry in Lemma |4.7 Left panel; 
cross-section at = T, with the initial rectangles rj (yellow) and fj (red) 
shown. Right panel: cross-section at = 0, with the hner structure of 
R 3 (-, 0 ) shown. 


such that each of them has aspect ratio no larger than 3, fj C rj C Ql, the rj are pairwise 
disjoint, \fj\ ~ pg, \rj \ ~ dg, with 


Assume 

(4.7) 


\fj\ G 27rN for all j. 


«Po > 1 


Then there is a pair admissible in the sense of Definition 2.5, and such that 

on Ql , 




(4.8) 


r^\x,B,QL,T\<Kp,TN^ + K 


2 , ,2Po^O^^ 


T 


and, with bj := K/\/2{\fj\/\rj\), 


(4.9) 


Ar2 


Bsi', 0 ) - ^ bjlr 

i=i 


< ^pIn^ + 


K^PgiV^do 


h-^/AQl) 

The set u := {x E Ql,t '■ x{^) = 1} is formed by the union of finitely many sheared 
parallelepipeds, with two faces normal to e^, and 


(4.10) 


\DX\ 

^ Jql.t 









































































































































For any z G (0, T) the flux of B across each connected component of {xa = z} is an 
integer multiple of 271. 


Proof. We first observe that a simple construction which obeys all kinematic constraints is 
obtained using a pattern which does not depend on X 3 . From the boundary data at X 3 = T 
we see that this necessarily is 


K 


B = ( 0 , 0 , 


X(x',X3) = Ir 




(X 


A simple computation shows that 

F'“[x,B,Ql.t] 


KpaTN^ 


and, using Lemma 4.4 


4(-,o) 


E 

t=l 


bilr 








At the same time |(c(;)i/k \ wI < N^T{pq/k, + 1/k^). If T < do, the proof is concluded. 
We observe that if upo < 1 the energy estimates would also hold, but not the one on the 
measure of In the following we assume T 3> do- 

We choose I eN such that 


(4.11) 


mm < Kpo 


h 7 

do 


Possibly reducing / by a few units, which does not affect the statement, we can assume 
that 


(4.12) 


4'< 


3 \ 

Stt 


K\r.i 


877^/2 


for all j 


(to see this, one observes that K\fj\ ~ Kpl > {Kpofl). The construction is performed 
independently in each set rj x (0, T), we take both helds to vanish outside the union of these 
sets. Let Rj,i,h be the rectangles given in Lemma 4.2 for 1 < ^ 1 < h < 4®, starting 


from rj, using the held i?* = bj. We denote by i the rehnement level, by h the numbering 
of the rectangles at each level. By construction we have \\Rj,i,h\ — 4“®|rj|| < jbj for all 


j,i, h] with (4.12) we obtain \Rj^i^h\ > 4 
and analogously \Rj,i,h\ < |4“ 


\rj\—A7r/bj=4: *|rj|(l —47r4®/(|rj|6j)) > 44 ®|r 




' 3 \- 


Therefore \R. 




We localize the magnetic held in an appropriate subset of each of the Rj,i,h- This is done 
by the family of inner rectangles {Rj^i^h)h=i,...p' which we now dehne: we let Rj^i^h have the 
same center and aspect ratio as and area given by 












Figure 4. Sketch of the vertical structure in Lemma [4.71 The thickness of 
the tubes at stage i is pi, their horizontal separation di, the distance between 
two vertical steps f*. 


We recall that bj is dehned so that bj\rj\ = K\rj\/y/2 G 27rM, in particular bj < K/y/2. 
Therefore, Rj^i^h C Rj,i,h, and the locally optimal held k/\/2 carries over Rj,i,h the same 
hux that bj carries over Rjj^h- At the level i = 0 this dehnition gives Rqj = fj. 

To estimate the size of the rectangles we dehne 


di := 2 *do and pi := po2 ®. 

Then \Rj,i,h\ ~ ~ d] and, correspondingly, \Rj,i,h\ = \Rj,i,h\ ~ Pi for all 

i = 1,... 

We now dehne the vertical structure, see Figure]^ The rehnement steps, labeled by i in 
the decomposition of the rectangles, will describe the structure at diherent levels, which 
are labeled j/j. Precisely, for some 6 G (0,1) chosen below, we dehne 


:= T0* 


U :=yi-y^+l = Te\l-^). 


1, is done using Lemma 4.6 


By (4.11), if 6* > 1/4 we obtain ti > di for all i < I. 

The explicit construction in x (pi+i, pi), for i = 0,..., /- 
At each step we interpolate between Rjj^h and the four corresponding Rjj^i^h^, Z = 1, 2, 3,4. 
Since the side length of Rjj^h is controlled by di, the side length of the inner rectangle where 
the held is located is controlled by pi, and di < ti, the energy in x (pj+i, yi) is bounded 
by 


Ar2 


j=l h=l 


W2' 

. 9 9 ^0 

KpiU + K — 


< N\poto{29y + 


pjdl 

toW ■ 














The series converges for all 6 G (1/4,1/2). Choosing 6 = 1/3 we conclnde that the energy 
in X {yi,T) is bonnded by 

I Ar 2 4 i 


df 


2=1 j=l h=l *- 


4. I 2 Z 2 

^Pi^i ^ Pi 

L‘ 


/ rr..r2 2 2 ^^ 0 ^^ 

< kpqTN + k Pq 


T 


We still have to consider the region Ql x (O,?//). Here we set 5' = 0 and (y, 0 : 3 ) = 
(x, B){x', Hi). The only energy contribntion comes from the snrface term, and is npiyj. 
Since yi < tj, this is controlled by the last summand in the series, and therefore does not 
change the scaling of the total energy. 


We hnally estimate the boundary term. Using Lemma [4.4| we obtain 

Ar2 4/ 


Ar2 4/ 

EE 

j=l h=l 


-^1 


72 


R. 






r < K- 


j=l h=l 


^3,1,h 


3/2 

~ 2 ^ 


We observe that 'Ylh^Rjih — ^y the construction of the rectangles. Inserting, from 
the definition (4.11), 2~^ ~ {kpo)~^ + d^/T, concludes the proof of (4.9). 

To prove (4.10) from Lemma 4.6 it suffices to show that 1 /k < p/, i.e., npi > 1. Since 


Kpi = 2 ^Kpo^ this follows immediately from (4.11). 


□ 


Proof of Theorem f.l. We £x fc G N, set N := 2^, and let {Rk,j}i<j<N^ be the rectangles 
given in Lemma 4.2 applied to i?o,i '■= Ql = (0, L)^ with := 6 ext- We assume 

(4.13) 

so that \Rk,j\ ~ Lf/N"^. We hrst localize the flux to the central part of each of the 
rectangles. This is the step in which the two regimes differ. Fix a factor 7 G (0,1], chosen 
below, and let be a rectangle with the same center and aspect ratio as Rkj, scaled by 
a factor 7 . We shall concentrate the entire flux over Rkj into the smaller rectangle rj, so 
that the magnetic held over rj is b := &ext/7^; from feextl-Rfcjl ^ 27rZ we obtain b\rj\ G 27rZ. 


Since 6 ext = J2j=l ^extlfifej , Lemma 


4.4 


yields 


(4.14) 


Ar 2 
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b^n ^extl 






rp/H 


7<1 




7<1 


The factor l-y<i represents the fact that this term is only present if 7 < 1. In the case 
7 = 1 , indeed, we have ^ = 1 on Ql, and therefore this term vanishes. 

with the given set of r^, do = 7 ^, k, L, N as above, and 


We then use Lemma 


4.7 


,_ ' j) /TV’ 

thickness T/2, on the set Ql x (0,T/2) (the other half is symmetric and not discussed 
explicitly). The inner rectangles are chosen so that fj has the same center and aspect ratio 
as rj and area given by K\fj\/\/2 = b\rj\, correspondingly the length scale is 

V 2 / 1/2 

K I \ K 


V2\ 


Po ■— do 


























Since fj needs to be a subset of r^, this is possible only if 6 < k/a/2, therefore it is only 
possible if 7 is chosen such that 


(4.15) 


bextV^'\ 


1/2 

< 7 < 1 • 


Since bext < ^/2, this set is non empty. At the same time, the condition (4.7) is satished 
provided that 


(4.16) 


(&extK)^/^— > 1 . 


To estimate the boundary term we combine (4.14) with (4.9) 

l|-S3(', 0) — 5ext|| u-i/ 




<2||i?3(-,0) btr-\\ h-^/2(o,) 

i=i 


+ 2 ||y^Wrj ^ext||/^-l/2(Q^) 
i=i 


and obtain 
(4.17) 


|i?3(-, 0) - 6ext||^-l/2(Q^) < PO 0 


T 


We then extend B to Ql x {— 00 , 0) by Lemma 4.5, and dehne y and B on Ql x (T/2, cx)) 
by symmetry, (y, i?)(a;', T/2 + z) = (y,T)(x',T/2 — z). Recalling (4.8) we see that the 
total energy is bounded by 

Fix. B] < .p,>TN^ + . 

Since po < the last term is smaller than the second one and therefore can be neglected. 
Dividing by the area and inserting the dehnitions po = Lb]I^/{N and do = L'y/N 
gives 


(4.18) 


TO ~ ^ ^ext T jY 27 ^ ’ 


T 2 


The construction is possible for all 77 > 1, 7 G (0,1] which obey (4.13), (4.15), and (4.16). 
To conclude the proof it suffices to choose these parameters appropriately. 


Choice of the parameters: intermediate regime. We assume here 

^ 5/7 


(4.19) 


"^ext > 


y2T2/7 


In this regime, we balance the hrst two terms in (4.18) by choosing the length scale L/N 


as 


W _ {KbextY^^ 

L ~ aT^C ’ 


7 = 1 


and 


N = inf{2^ : fc G N, 2^ > W} where 




















where a is a number of order 1 chosen below. We assume W > 1 so that the rounding of 
W to the next power of 2 does not modify it by more than a factor of 2, N^: < N < 2W- 
One obtains 


F[X,B] 

L 2 




where we used that one term disappears because 7 = 1 ; the 6 ext term can be dropped since 

= (6ext/«:)'/'(«T)-V3 < 1. 

It remains to check that th e choi ces made are admissible. Condition (4.16) translates 
into a{KbextT‘^YF > 2 . Using (4.19), kT > 1 and assuming a > 4 one can easily see that 
it is satisfied. Since N < 2W, condition (4.13) translates into 


^2 2^4 

^ext~^ 

K 


1/3 


327r 


> 

a2 


Using hrst (4.19) and then kT > 1 and k < 1/2 one can easily check that the parenthesis 
is at least 4, therefore it suffices to choose a = 8 . Condition (4.15) is immediate. Finally, 
we check that > 1. This is equivalent to 


L > 


^2/3 


{nbexty/^^ • 


Choice of the parameters: extreme regime. In this case we assume 

^5/7 


(4.20) 


^ext ^ 


v/2T2/7 ■ 


In this regime, we can actually make the three first terms in (4.18) balance by choosing 
the length scale L/iVand 7 according to 


7 


= 2 V 4 /LJ^ and 


K 


5/14 


N = inf{ 2 ^ : fc e N, 2 ^ > W} 


7/2 


-^ext 


where — = ■ , „, , 

L ’ 


for some a > 0 chosen below. Again, we require > 1 so that < N < 2N^. This 
gives do ~ and 

since kT > 1. 


We turn to checking that the choices made are admissible. The assumption (4.20) is 
equivalent to 7 < 1. The other inequality in (4.15) is fulfilled by kT > 1 (this is the r eason 
for inserting the factors 2 }^ and a/2 in the dehnition of 7 and ( 4.20 )). Condition (4.16) 
becomes airCryC > 2 , which is true for any a > 2. Condition (4.13) becomes 


c? [K 


TfC 


> 1 


3277 K 
























which again is satished if a = 8 . The fact that > 1 translates into 


L > 


a/2 

"'ext 


This conclndes the proof. 


□ 


4.2. Construction for the Ginzburg-Landau functional. We hnally give the npper 
bonnd construction for the Ginzburg-Landau functional. We start from the constructions 


given in Section 4.1 for the sharp-interface functional F. 


The hrst step is to construct the vector potential A from the magnetic held B. We use 
the following lemma, which is a variant of Hodge’s decomposition in the current geometry. 


Lemma 4.8. Let B E 

^'qlx'R 

(4.21) 


loc 


Ql periodic, such that div B = 0 distributionally and 


\B\‘^dx < oo. Then there is A & also Qi-periodic, such that 


Proof. From Lemma 2.2 we obtain 


V X A = B. 


Hs dx' = 0 for all G 


Jq{z) 

(in the sense of traces). We dehne 


B{x) := B{x) - — 


Bdx' 




SO that all components of B have average zero on each cross section Q{z). Notice that 
.03 = 03 , and V'0 = V'0. This implies 0 G 0^, and div0 = 0. 

The rest of the proof gives a construction of A such that B = V x A. If one could solve 
directly AT = 0 with divT = 0, then V x T would work. Since we are working in an 
unbounded domain with mixed boundary conditions, for completeness we give an explicit 
construction of A based on Fourier series. 

We Fourier transform to obtain coefficients b{k) such that 


B(x) = 


j dks 

k'£2nZ^/L 


The transformation we just performed ensures that b{k) = 0 whenever k' = 0 (this is the 
reason to consider 0 instead of 0). Further, k ■ b{k) = 0 for all k. 

We dehne 

{ —ik X h{k) 

P 

0 if fc' = 0 , 

so that for all k 


if fc' ^ 0 


ik X d{k) = h{k). 






The family a{k) also corresponds to a converging Fourier series, since \a{k)\ < L\b{k)\/{2n) 
(this is the step where it is important that we requested k' ^ 0 , and not merely k ^ 0 ). 

In real space, we set 

A{x) := f dk 3 • 

fc'e27rZ2/L 

Clearly A G (M^; M^), it is Q^-periodic and obeys 

V X A = B, 

both sides being in L‘^{Ql x We hnally set 


A{x) = A{x) - 


(cs X B)dx. 


' QlX (0,X3) 

The correction depends only on X 3 . Therefore, recalling -Bsdx' = 0, 

V X A = V X A — — I 63 X {63 X B)dx 

^ JQ(x3) 

= B + — I B dx' = B . 

B Jq(x3) 

This concludes the proof. 

Theorem 4.9. For any 6ext, k, L,T > 0 such that 

1 

and 


□ 


4&ext 


kT > 1, 


L sufficiently large (in the sense of (4-2) and be^tL"^ G 27rZ there is a pair {u,A) G Tfper, 
E[u, A] < min 


such that 


Proof. Let x and B be the functions constructed in Theorem |4.1[ Starting from B, we 
obtain A from Lemma 4.8 Consider now the superconducting domain ca := {x G Ql,t '■ 
x{x) = 0}. Here H is a curl-free vector held, hence it is locally the gradient of some 
potential 6 . The domain is multiply connected, but the hux of B across each tube is an 
integer multiple of 27r, hence we can globally write A as the gradient of a multi-valued 
function 6 , such that 6 mod 27r is single-valued. 

We set p = 0 in the normal phase, and let it grow to 1 in the superconducting phase, on 
a length scale 1 /k, 

p{x) ;= min{l, dist^(a;, ca)}. 

The distance function is understood, as usual, Q^-periodic in the hrst two components. 
Finally, we set 

u{x) := . 

Since V 6 = A whenever p 7 ^ 0, we have 

\Vau? = |Vp'/T- 





The B' part of the energy is identical, and so is the onter held. It remains to treat the 
conpling term. In u we have p = 0 and x = 1) hence 

Outside oj we have i? = 0, hence 


The hrst term is exactly the one appearing in F. Therefore the E[u,A\ < F[x,S] + Es, 
where 

Es := [ + K^{x-{1 - p)y] dx . 

d Ql,t 

Let {u)i/k be a 1 /k- neighbourhood of u. Then in uj we have y = 1 = 1 — p, outside {uj)i/k 
we have y = 0 = 1 — p, and recalling |Vp^/^| < k we obtain 


Es < 2 k^\{uj)i/^\uj\ . 

Recalling (4.3) we conclude Es ^ n Iqlt This concludes the proof. 


□ 
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